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1. Introduction 

Experiments at the Large Hadron Collider (LHC) will allow to test standard model at very 
high energies. Here we are interested in Quantum Chromodynamics (QCD) processes like forward 
jet production and interplay of QCD with Electro-Weak interactions which lead to W, Z bosons 
production in hadronic collision. The first of listed processes is of interest since it will allow for 
better understanding of partonic structure of the proton at extreme energies where new phenomena 
like saturation |Jlj, £| [||] may occur. The other process is an important candle of Electro-Weak the- 
ory which has to be tested at various energy scales in order to have full understanding of it. This 
process will also serve as one to calibrate calorimeters and to measure luminosities. 
The large center of mass energy at the LHC will require a need to apply QCD resummation ap- 
proaches capable to account for multiple scales in the problem. Namely, one has to account for 
logarithms of type a $ ln(jU / Aqcd) where /I might be pr of jet or mass of W/Z and logarithms 
coming from the fact that at least parton densities in one of the incoming proton will be probed 
at very small longitudinal momentum fraction x giving rise to logarithms of type a"ln m 1/x The 
theoretical framework to resume consistently both kinds of logarithmic corrections in QCD calcu- 
lations is based on high-energy factorization at fixed transverse momentum. This formulation 
depends on unintegrated distributions for parton splitting, obeying appropriate evolution equations, 
and short-distance, process-dependent matrix elements. The unintegrated-level evolution is given 
by evolution equations in rapidity, or angle, parameters. Different forms of the evolution, valid 
in different kinematic regions, are available, see [|5[ §J, and references therein, for recent work in 
this area and reviews. In this article we present recently obtained results for hard matrix elements 
needed in application of factorisation formulas for both of introduced above processes. In Sec. 2 
we introduce the basic structure of jet production in the LHC forward region. In Sec. 2.1 we con- 
sider associated parton showering effects. In Sec. 2.3 we consider effects from the short-distance 
matrix elements that control the resummation of logarithmically enhanced corrections in yfs/Ej, 
where Ej is the hard jet transverse energy. In Sec. 3 we briefly discuss central production of W, 
Z bosons focusing on observables which are sensitive to different showering methods and different 
assumptions on initial state of colliding partons. We give concluding remarks in Sec. 4. 



2. Forward jets 

The hadroproduction of a forward jet associated with hard final state X is pictured in Fig. [I]. 
The kinematics of the process is characterized by the large ratio of sub-energies S2A1 3> 1 and 
highly asymmetric longitudinal momenta in the partonic initial state, q& ■ ps 3> qB • Pa- At the 
LHC the use of forward calorimeters allows one to measure events where jet transverse momenta 
Px > 20 GeV are produced several units of rapidity apart, Ay > 4 4- 6 [{7|, ]1| Q]. Working at polar 
angles that are small but sufficiently far from the beam axis not to be affected by beam remnants, 
one measures azimuthal plane correlations between high-p^ events (Fig. ^ widely separated in 
rapidity [g 0. 

The presence of multiple large-momentum scales implies that, as recognized in [11, |12[ p"3| ], 
reliable theoretical predictions for forward jets can only be obtained after summing logarithmic 
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Figure 1: Jet production in the forward rapidity region in hadron-hadron collisions. 



QCD corrections at high energy to all orders in a s . This motivates efforts [ ]14] , \15\ |lq , |17[ ] to 
construct new, improved algorithms for Monte Carlo event generators capable of describing jet 
production beyond the central rapidity region. 




Figure 2: (Left) High-px events in the forward and central detectors; (right) azimuthal plane segmentation. 

In the LHC forward kinematics, realistic phenomenology of hadronic jet final states requires 
taking account of both logarithms of the large rapidity interval (of high-energy type) and loga- 
rithms of the hard transverse momentum (of collinear type). The theoretical framework to resume 
consistently both kinds of logarithmic corrections in QCD calculations is based on high-energy 
factorization at fixed transverse momentum [fill. 



Ref. [ |10| ] investigates forward jets in this framework. It presents the short-distance matrix 
elements needed to evaluate the factorization formula, including all partonic channels, in a fully 
exclusive form. On one hand, once convoluted with the BFKL off-shell gluon Green's function 
according to the method of [ft], these matrix elements control the summation of high-energy loga- 
rithmic corrections to the jet cross sections. They contain contributions both to the next-to-leading- 



order BFKL kernel [18] and to the jet impact factors [19, 20]. On the other hand, they can be used 



in a shower Monte Carlo generator implementing parton-branching kernels at unintegrated level 



(see e.g. [21 , 22] for recent works) to generate fully exclusive events. 

The high-energy factorized form [fj, [T(| |l9|] of the forward-jet cross section is represented in 
Fig. ^]a. Initial-state parton configurations contributing to forward production are asymmetric, with 



'Analogous observation applies to forward jets associated to deeply inelastic scattering [ p3[ ^j ]. Indeed, measure- 
ments of forward jet cross sections at Hera J2^ ] have illustrated that either fixed-order next-to-leading calculations or 
standard shower Monte Carlos p^ , [l4j], e.g. PYTHIA or HERWIG, are not able to describe forward jet ep data. 



3 



Studies of forward jets and production ofW,Z 



the partem in the top subgraph being probed near the mass shell and large x, while the partem in 
the bottom subgraph is off-shell and small-x. The jet cross section differential in the final-state 
transverse momentum Q± and azimuthal angle (p is given schematically by [Q, [l9|] 

where (8) specifies a convolution in both longitudinal and transverse momenta, a is the hard scatter- 
ing cross section, calculable from a suitable off-shell continuation of perturbative matrix elements, 
ty a / A is the distribution of parton a in hadron A obtained from near-collinear shower evolution, and 
ty g * i B is the gluon unintegrated distribution in hadron B obtained from non-collinear, transverse 
momentum dependent shower evolution. 




(a) (b) 



Figure 3: (a) Factorized structure of the cross section; (b) a typical contribution to the qg channel matrix 
element. 

In the next section we comment on the initial-state shower evolution. In Sec. 2.2 we turn to 
hard-scattering contributions. 

2.1 Parton shower evolution 

Parton distributions can be obtained by parton-shower Monte Carlo methods via branching al- 



gorithms based on collinear evolution of the jets developing from the hard event [|27p. The branch- 
ing probability can be given in terms of two basic quantities (Fig. ^|), the splitting functions at the 
vertices of the parton cascade and the form factors to go from one vertex to the other. An important 



ingredient of this approach is the inclusion of soft-gluon coherence effects Q27| , |28| , [29| ] through 
angular ordering of the emissions in the shower. 

Corrections to collinear-ordered showers, however, arise in high-energy processes with multi- 
ple hard scales [|9, 3C, 31], as is the case with the production of jets at forward rapidities in Fig. [[]. 
In particular, new color-coherence effects set in in this regime due to emissions from internal lines 
in the branching decay chain [Q, [l^, [32|] that involve space-like partons carrying small longitudi- 
nal momentum fractions. The picture of the coherent branching is modified in this case because 
the emission currents become dependent on the total transverse momentum transmitted down the 
initial-state parton decay chain 



19, 30, 31 



]. Correspondingly, one needs to work at the level 
of unintegrated splitting functions and partonic distributions ||5|, 0] in order to take into account 
color coherence not only for large x but also for small x in the angular region (Fig. ^J) 



a/x > oci > a 



(2.2) 
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Figure 4: Parton branching in terms of splitting probabilities and form factors. 

where the angles a for the partons radiated from the initial-state shower are taken with respect to 
the initial beam jet direction, and increase with increasing off-shellness. 



P 

+ vl + 



Figure 5: (left) Coherent radiation in the space-like parton shower for x -C 1 ; (right) the unintegrated split- 
ting function including small-x virtual corrections. 

The case of LHC forward jet production is a multiple-scale problem where coherence effects 
of the kind above enter, in the factorization formula dOh, both the short-distance factor a and 



the long-distance factor 0. Contributions from the coherence region ( |2.2| ) are potentially enhanced 
by terms q"ln m \/s/px_ where y/s is the total center-of-mass energy and p± is the jet transverse 
momentum 2 . These contributions represent corrections to the angular ordering implemented in 
collinear showers and are not included at present in standard Monte Carlo generators [|27|]. Work 
to develop methods for unintegrated shower evolution, capable of including such corrections, is 
underway by several authors. 



The proposal pip incorporates NLO corrections to flavor non-singlet QCD evolution in an 
unintegrated-level Monte Carlo. The approach is based on the generalized ladder expansion of [[35|], 
which is extended to the high-energy region in [p6Q. This approach could in principle be applied 
generally, including flavor singlet evolution, and used to treat also forward hard processes. 

2.2 The factorizing hard cross sections 

Logarithmic corrections for large rapidity y ~ In s/p\ are resummed to all orders in a s via 
Eq. (2A), by convoluting (Fig. |3j> unintegrated distribution functions with well-prescribed short- 



2 Terms with m > n are known to drop out from inclusive processes due to strong cancellations associated with 
coherence, so that, for instance, the anomalous dimensions f J for space-like evolution receive at most single-logarithmic 
corrections at high energy jll], pq ]. This need not be the case for exclusive jet distributions, where such cancellations 
are not present and one may expect larger enhancements. 
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distance matrix elements, obtained from the high-energy limit of higher-order scattering ampli- 
tudes flio[ yjl. With reference to Fig. ||b, in the forward production region we have (/?4 + p(,) 2 3> 
(P3 + Pa) 2 and longitudinal momentum ordering, so that 



Pe 



(2.3) 



Here E,\ and £2 are longitudinal momentum fractions, and is the di-jet transverse momentum in 
the laboratory frame. It is convenient to define the rapidity-weighted average Q± = (1 — v)p±4 — 
Vp±$, with v = (j)2 ■ P<x)/p2 • (pi — P5). In Fig. ||b Eq. (7A) factorizes the high-energy qg amplitude 
in front of the (unintegrated) distribution from the splitting in the bottom subgraph. The factoriza- 
tion in terms of this parton splitting distribution is valid at large y not only in the collinear region 
but also in the large-angle emission region As a result the rapidity resummation is carried out 
consistently with perturbative high-i2_L corrections [Q, [1^] at any fixed order in a s . 
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Figure 6: The kj/Qr dependence of the factorizing qg hard cross section at high energy Jl0|]: (left) C] 
term; (right) CfCa term. 



The explicit expressions for the relevant high-energy amplitudes are given in p0|]. Figs, g and 
[7] illustrate features of the factorizing matrix elements, partially integrated over final states. We plot 
distributions differential in Q± and azimuthal angle (p (coscp = Q± • k±/\Q±\ \k±\) for the case of 
the qg channel. Fig. |6] shows the dependence on k±, which measures the distribution of the third jet 
recoiling against the leading di-jet system. Fig. [7] shows the energy dependence. 

C A C F term qg channel 

& T dfr 



C F term qg channel 



5x10 



3x10 



1x10 




4x10 



3x10 



2x10 



1x10 




9 10 s ^ 



5 6 7 8 9 10° 
Figure 7: The energy dependence of the qg hard cross section [ |To| ] (kj /Qt = 1). 

The region k±/Q± — > in Fig. |6| corresponds to the leading-order process with two back- 
to-back jets. The resummation of the higher-order logarithmic corrections for large y ~ lns/p 2 ± 
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is precisely determined [Q, [19|] by integrating the u-pdfs over the ^-distribution in Fig. ^. So 
the results in Fig. ^ illustrate quantitatively the significance of contributions with k±_ ~ <2_l in the 
large-j region. The role of coherence from multi-gluon emission is to set the dynamical cut-off 
at values of k± of order Q±. Non-negligible effects arise at high energy from the finite-^ tail. 
These effects are not included in collinear-branching generators (and only partially in fixed-order 
perturbative calculations), and become more and more important as the jets are observed at large 
rapidity separations. The dependence on the azimuthal angle in Figs. ^ and ^ is also relevant, as 
forward jet measurements will rely on azimuthal plane correlations between jets far apart in rapidity 
(Fig.0). 

Results for all other partonic channels are given in [[K]]. After including partem showering [Bq], 
quark- and gluon-initiated contributions are of comparable size in the LHC forward kinematics: re- 
alistic phenomenology requires including all channels. Note also that since the forward kinematics 
selects asymmetric partem momentum fractions, effects due to the x — ► 1 endpoint behavior [Q] at 
the fully unintegrated level may become relevant as well. 



3. W, Z boson production 

Central production of W/Z bosons at LHC energies will be dominated by gluonic compo- 
nent of partonic system. The longitudinal components of incoming gluons four momenta will 
be small and of the same order what will result in the central production. As in the forward jet 
case the framework to consistently account for this kinematical set up is provided by k± factori- 
sation approach where o(g*g* — > q(W /Z)q) is an off shell continuation of hard matrix element 



for o(gg — ► q(W/Z)q) [ [39 , 40] which allows to include effects coming from finite transversal mo- 
menta of gluons. The interesting observable to calculate for this production process is the angular 
distance between highest p± jet and p± of Z or W, which allows for estimation of uncertainties 
of theoretical predictions. On Fig. |9] (left) we show comparison of calculation based on k± fac- 



torisation approach Q39| , |40p and collinear approach at LO and NLO. We see that the distributions 
are considerably different for small angles in case of k± factorisation approach and collinear LO. 
While NLO collinear is quite similar to the one obtained in k± factorisation approach. The reason 
for that comes from the momentum conservation in LO collinear approach which forbids events in 
region from to n/2. This not the case in k± factorisation approach and NLO collinear calculation 
where the additional transversal momentum flow allows momenta of Z, b and b to be unbalanced. 
The k± factorisation formula while applied to this process takes form: 

a = fo /A ®d(g*g* ^q(W/Z)q)<j> b y B , (3.1) 

The other important distribution is the differential cross section in the transversal momentum of Z 
bb system || (right). This distributions are rather different at small p zh j,^ region. The difference 
is due to the fact that in the k± factorisation approach subleading terms of all orders from point of 
view of collinear approach are taken into account. We also see that at higher values of p zfc £i both 
approaches agree. 
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Figure 8: Jet production in the forward rapidity region in hadron-hadron collisions. 







10' 


[ 






cofTmear LO 






colTmear LO 


1 


collinear NLO 






collinear NLO 




: — kt 


e 




— kt 






° 3 






S3 




n 10" 3 
c 






J 5 

< 




If . 






■o 

o 10 
■o 




Q~10- 

■o 

"6 










~° 10 s 










10 s 

















0.5 1 1.5 2 2.5 3 50 100 150 200 250 300 350 400 

A *Zhb P zb Bi 

Figure 9: Distributions of the distance in azimuthal angle of Z and highest px quark or antiquark. Calcula- 
tion with massive b-quarks (left). Comparison of cross sections in transverse momentum of the produced Z 
gauge boson. Calculation with massless ^-quarks (right) 

4. Conclusions 

Forward + central detectors at the LHC allow jet correlations to be measured across rapidity 
intervals of several units, Ay > 4 -j- 6. Such multi-jet states can be relevant to new particle dis- 
covery processes as well as new aspects of standard model physics. Existing sets of forward-jet 
data in ep collisions, much more limited than the potential LHC yield, indicate that neither conven- 
tional parton-showering Monte Carlos nor next-to-leading-order QCD calculations are capable of 
describing forward jet phenomenology. Improved methods to evaluate QCD predictions are needed 
to treat the multi-scale region implied by the forward kinematics. In this article we have discussed 
ongoing progress, examining in particular factorization properties of multi-parton matrix elements 
in the forward region, and prospects to include parton-showering effects with gluon coherence not 
only in the collinear region but also in the large-angle emission region. We also have discussed 
predictions for W and Z central production based on different showering schemes. 
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